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Abstract. One of the main sources of gravitational waves for the LISA space- 
borne interferometer are galactic binary systems. The waveforms for these sources 
are represented by eight parameters of which four are intrinsic and four are 
extrinsic to the system. Geometrically, these signals exist in an 8-d parameter 
space. By calculating the metric tensor on this space, we calculate the number of 
templates needed to search for such sources. We show in this study that below a 
particular monochromatic frequency of /o ~ 1.6 X 10~^ Hz we can ignore one of 
the intrinsic parameters and search over a 7-d space. Beyond this frequency, we 
have a change in dimensionality of the parameter space from 7 to 8 dimensions. 
This sudden change in dimensionality results in a change in the scaling of template 
number as a function of the monochromatic frequency from ~ /g'^^ to ~ /g **. 
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1. Introduction 

One of the main sources of gravitational waves for the planned space-borne Laser 
Interferometer Space Antenna (LISA) are quasi-monochromatic galactic binary 
systems. The LISA detector is composed of three spacecraft which form an equatorial 
triangle and will work in the frequency range 10^^ < f /Hz < 1. The center of mass 
of the constellation traces out a circular orbit around the Sun at a distance of 1 AU 
and lies about 20° behind the Earth. The three spacecraft cartwheel in retrograde 
motion as they move around the Sun with a period of one year. This motion induces 
amplitude, frequency and phase modulations in the gravitational wave signalUE]- The 
amplitude modulation is caused by the antenna pattern being moved across the sky. 
With the three arms, LISA acts like a pair of two-arm detectors. The constellation is 
set up such that the arm lengths of the detector are equal with an angle of 7r/3 between 
each detector arm. This is equivalent to the response of two 90° interferometers rotated 
by an angle of 7r/4 radians with respect to each other. Because LISA can be thought 
of as two separate detectors, measuring different polarizations of the GW, the phase 
modulations are cause by combinations of the polarizations. Finally, the frequency 
or Doppler modulations are caused by the motion of the detector with respect to the 
source. The noise in the LISA detector is complicated by the inclusion of a transfer 
function at high frequencies. Due to this, we will work, in this study, in the Low 
Frequency Approximation (LFA). In the LFA we can set the transfer function due to 
the arm length of the detector equal to unity for / <C /* = 10~^ Hz. 

In order to detect galactic binaries one can employ the method of matched filtering 
P]. Briefly, the method works as follows: Firstly, one creates a set of waveforms, or 
templates as they are called, that depend on a number of parameters of the source 
and its location and orientation relative to the detector. These templates are then 
cross-correlated with the detector output weighted by the inverse of the noise spectral 
density. If a signal, whose parameters are close to one of the template waveforms, 
is actually present in the detector output then the cross-correlation builds up, with 
the dominant contribution coming from frequencies where the noise spectral density is 
low. Thus, in the presence of a sufficiently strong signal the correlation will be much 
larger than the RMS correlation in the absence of any signal. How large should it be 
before we can be confident about the presence of a signal depends on the combination 
of the number of sources and the false alarm probability. The effectiveness of matched 
filtering depends on how well the phase evolution of the waveform is known. Even 
tiny instantaneous differences, as low as one part in 10"^ in the phase of the true signal 
that might be present in the detector output and the template that is used to dig it 
out could lead to a cumulative difference of several radians since one integrates over 
several hundreds to several thousands of cycles. 

Using a geometric method U, we can treat the parameters of the systems as a 
set of basis coordinates. We can then define a scalar product and metric tensor in this 
vector space. Having done this we can show that we only need to carry out our search 
over a lower dimensional sub-space. The search for galactic binaries will then involve 
the placement of a grid of templates in this sub-space. 

2. The Gravitational Waveform. 

The strain of the gravitational wave with both polarizations is 

h{t) ^ h+ F+ + hy^ , (1) 
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where the two polarizations are given by 

h+ (t) = ^0 (1 + cos^ i) cos ($ (t) - (fio) , 
hx {t) — —2 ^0 cos L sin ($ (t) — i^o) , 



(2) 
(3) 



where Aq is a constant initial amplitude, t the angle of inclination of the source, and 
(fa is a constant initial phase. The phase of the gravitational wave, for a circular 
orbit is defined by 

- 27r/o< + 7r/o<2 + 27r (/o + /oi) R sin(0) cos (2^/„< - , (4) 
where fo is the quasi-monochromatic frequency and /q is its first time derivative given 

by 



^ = ^7r^/3_A/(5/3jn/3^ 
dt 5 



(5) 



Here, A4c = mrf'/^ ^ is the chirp mass, and, m = mi + m2 and rj = mim2/m^ are the 
total mass and reduced mass ratio. (6*, </>) define the location of the source, /,„ — 1/T 
is the modulation frequency and T is the number of seconds in a year, and R is the 
light travel time for 1 AU 500 sees). The quantities and are defined in the 
LFA by IHl 

F+it; e, A) = i [cos(2V')i^+(i; 9, 0, A) - sin(2V^)i?>< {t- 6, ^, A)] , (6) 

F"" (t; V, 61, 0, A) = i [sin(2V')£'+(t; 61, 0, A) + cos(2V')i:''' (t; 6*, 0, A)] , (7) 

where ij) is the polarization angle of the wave and A = or 37r/2 defines the two-arm 
combination of LISA from which the strain is coming. The detector pattern functions 
are given by 



D+{t) = — - aesin^r^) sin(2a(<) - 2A) + (3 + cos(2e')) 
64 L 

cos(2(?!)) I 9 sin(2A) - sin(4Q:(t) - 2A) | -1- sin(2(/)) | cos (4a(t) - 2A) 9 cos(2A) | 
-4V3sin(26') ( sin(3a(i) -2\-(j)) - 3sin(a(t) - 2A + ( 



(8) 



D^'U) = 4: [ V3cos(6l) f 9 cos(2A - 2</>) - cos(4am - 2A - 20)) (9) 
16 L \ / 

-6 sin(0) ( cos(3a(t) - 2A - </>) -I- 3 cos(Q!(t) - 2A -I- 0)) , 
where a{t) ~ 2Tit/T is the orbital phase of the center of mass of the constellation. 



3. The geometric method. 

The set of N parameters which define the system can also be thought of as defining 
an M dimensional manifold. This template manifold shares the properties of a Hilbert 
space. On this manifold is it possible to define a scalar product of the two time series 
waveforms h{t) and g{t) as 

= p -^h{f)r{f) + c.c., (10) 
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where the tilde denotes the Fourier transform h{f ) = h{t) exp{2'Kift)dt of the 
time series waveform h{t), a * denotes complex conjugate and the c.c. denotes a 
complex conjugate term. The quantity Sn{f) represents the two-sided noise spectral 
density of the detector, which is given by [5] 

SM) = (/) + 8 (^1 + cos^ (^-^)) ^^^7]?' (11) 

where S^^{f) and S'^^'^'^^{f) are the position and acceleration noise respectively. As 
we are searching for monochromatic binaries, we can greatly simplify the calculation 
of the inner product. In this case, the noise at the quasi-monochromatic frequency 
can be treated as a constant and taken outside the integral. Then, due to Parseval's 
theorem, we can calculate the scalar product in the time domain. This simplifies 
Equation (fTn|l to 

m ^ f dth{t)g{t), (12) 

JnUO) Jo 

where T is the period of integration (which in this study we will assume to be T = 1 
year). If we consider unnormalized waveforms, we can write the metric tensor on this 
manifold in the form 

_ {d^h\d,h) {h\d^h) {h\d,h) 

{h\h) ih\h)' ■ (^'^ 

The first term on the right hand side is recognized as being the Fisher information 
matrix. The extra term on the right hand side arises from the fact that we have 
not pre-normalized the waveforms. The number of templates needed to carry out 
a gravitational wave search is simply the proper volume of the parameter space, 
divided by the proper volume of the template. Assuming a D-dimensional hyper- 
cubic template, the number of templates is given by 

^-hi^^ (14) 

where g = de^lg^^l is the determinant of the metric tensor, A is the parameter set 
which defines the manifold and dl^ is the proper volume of the template. While the 
galactic binaries are defined by an 8-d parameter set, we only need to search over 
four intrinsic parameters. The reason for this will be explained in the next section. 
In order to do this, we project from the initial 8-d space to the 4-d search sub-space 
using the recursive projection 

7y = grj , (15) 

gkk 

which carries out a projection onto the sub-space orthogonal to the fc-dimcnsion. It can 
be shown that this projection onto orthogonal sub-spaces is equivalent to maximizing 
over the projected parameter [7]. Using the fact that the proper distance between two 
templates is related to a pre-determined threshold, called a minimum match, MM, via 

MM 1 ~ ids^ = 1 - ig^^dA^dA'', (16) 
we can now write the proper volume of the template in D dimensional space as 



(17) 
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4. Results. 

For numerical purposes it simplifies things if we use the parameters p = /qT and 
q{p) ~ foT^- For quasi-monochromatic binaries, the initial parameter space is 
defined by the set A'^ = {Ao,i,(fQ,'il;,9,(j),p,q}. The first four of these parameters 
are kinematic or extrinsic parameters. The four parameters {9, (f>,p, q} are dynamical 
parameters dependent on the astrophysical system in question and are called intrinsic 
parameters (the parameters {9, (f>) are considered intrinsic due to the motion of the 
detector with respect to the source). 

Instead of having to search over the full 8-d search space, we only need to search 
of the sub-set of intrinsic parameters = {9,(t),p,q}. This is due to the fact that 
we can use the F-statistic to find the optimal values of the four extrinsic parameters. 
This is done as follows. We first write the strain of the gravitational wave in the form 

4 

h{t) = a, {Ao, I, ifio, V-) A' {t; 9, (j),p, q) . (18) 

i=l 

We then define four constants A^' (s |^') , where s is the signal, and find a solution 
for the Oi's in the form 

a, = (A/-i)^^. 7V^ (19) 

where 

M'^ = {A' \A^) . (20) 
If we now substitute the above into the expression for the reduced log-likelihood, i.e. 

In c {xn = {s \h (xn) - i {h {xn \h (xn) , (21) 

we obtain the F-statistic |H] 

.F=i(M-i)^^.7VW, (22) 

which automatically maximizes the log-likelihood over the extrinsic parameters and 
reduces the search space to the sub-space of intrinsic parameters. 

In calculating the number of templates needed, we take different limits of the 
parameter p, we set a maximum value of q corresponding to a maximum chirp mass 
of A4c = 20 Mq (giving us a distribution of sample values of q for each value of p) 
and use these as a bounding box for a Monte Carlo integration to calculate the proper 
volume of the search manifold. As the boundary of our search space is quite regular, 
we can approximate the proper volume of the search space as 

d^X^g^^^^^, (23) 

where V = A9A(f)ApAq is the volume of the search space, ^/gl is the value of the 
metric determinant at a particular point, and N is the total number of points used. 
We can see from Figure that choosing a maximum value of q{p) for each p provides 
us with a range of values to sample over during the Monte Carlo integration. At the 
same time as carrying out the Monte Carlo integration, we also compare the proper 
distance in the g-dimension for each value of p, i.e. 

Sq^ dq ^/g^, (24) 
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Figure 1. A log-linear plot of p against q{p) for galactic binaries. During the 
Monte Carlo integration we sample from the shaded area, where the boundary 
of this area corresponds to a maximum chirp mass of 20 Mq. The dashed line 
corresponds to the value of p beyond which we have to increase the dimensionality 
of the search space and consider the parameter q. 

with the proper distance between two templates in the g-dimension, dl, as defined 
by Equation 1171 It is only when Sq > dl that we need to take this parameter into 
account and increase the dimensionality of the initial search manifold from 7-d to 8-d. 
Increasing the dimensionality of the search manifold infers a corresponding increase in 
the dimensionality of the sub-manifold of intrinsic parameters from 3-d to 4-d. This 
increase happens at p* ^ 5 x 10* and we now have Sq ^ dl. This transition point is 
indicated by the dashed line in Figure (Q. 

In Figure (0) we have plotted the cumulative number of templates, J\f{fo), needed 
to carry out a matched filtering search for individual quasi-monochromatic galactic 
binaries up to a particular choice of /q. We have assumed a minimal match in this case 
of MM — 0.97. We can see from the plot that at a frequency of Jq — ^ 1.6x 10~^ 
Hz, there is a point of inflection and a sudden increase in the growth of template 
number. It is because of the increase in dimensionality that we have the sudden 
increase in the growth of template number. We should also point out that the position 
of the change of dimensionality is dependent on our choice of MM. A higher choice 
of MM will push the point corresponding to the change in dimensionality to the left 
on our graph. 

In terms of how the number of templates scales against frequency, this increase 
in dimensionality causes a big change. Below /q, the number of templates scales as 
~ /o'^^. After the change in dimensionality, this changes to ~ fo'^^- 



Figure 2. The cumulative number of templates needed to search for individual 
quasi-monochromatic galactic binaries up to a particular value of fo with the 
LISA interferometer. We assume in this figure that the observation period T = 1 
year. 



5. Conclusions 

In this study we liave calculated the number of templates needed to carry out 
the search for individual quasi-monochromatic galactic binaries with the LISA 
interferometer. Using the standard geometric method, we have treated the waveform 
parameters as the basis coordinates in an 8-d space. As most of the parameters are 
extrinsic to the system, we only need to carry out the search on a 4-d sub-manifold. 
We find that below a frequency of /g = 1.6 x 10"'^ Hz, the proper distance between 
two templates in the fo direction is greater than the proper coordinate distance in 
that direction. As a consequence, below this frequency we can completely ignore the 
parameter defining this direction and our initial search space reduces to 7-d. It is only 
above /q that the proper coordinate distance in that direction becomes equal to or 
greater than the proper distance between two templates. We then have an increase in 
the dimensionality of the search space from 7-d to 8-d. It is because of this increase 
in dimensionality that we have the increase in the scaling law from ~ /o '^^ to ~ fo '^^- 
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